Dissipation of radial oscillations in compact stars 
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We demonstrate that there exists a new mechanism for dissipating the energy of stellar oscillations. 
For neutron stars, in particular, we show that the mechanical energy of density perturbations is not 
only dissipated to heat via bulk viscosity, but also that energy is radiated away via neutrinos. This 
energy dissipation will be associated with a viscosity coefficient, the radiative viscosity, which is 
larger than the bulk viscosity in the case of non-strange quark matter and nuclear matter. 
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Introduction — Dissipative processes play an important 
role in the evolution of neutron stars. These processes are 
governed by transport coefficients, such as the heat and 
electrical conductivities, and the shear and bulk viscosi- 
ties. The conductivities are important for stellar cooling 
as well as for the magnetic field decay. Shear viscosity 
damps differential rotation in a star and, thus, leads to a 
uniform rigid-body rotation. Bulk viscosity, on the other 
hand, damps the density oscillations of the stellar mat- 
ter. Both differential rotation and oscillations could be 
excited in newly formed (hot) neutron stars, or could de- 
velop in old (cold) stars due to external perturbations, 
e.g., such as matter accreted from a companion star. 

The viscosity of nuclear and mixed phases of dense 
baryonic matter has been calculated under various condi- 
tions and assumptions over the last three decades [l|-[lot . 
The bulk viscosity of normal conducting strange quark 
matter was also calculated The latter might be 

relevant if the baryon density in the central regions of 
neutron stars is so high that matter becomes deconfined. 

The main goal of this letter is to re-investigate the 
dissipation of density oscillations of high-density matter 
present in the interior of neutron stars. Using two-flavor 
quark matter as an example, it will be shown that the 
density perturbations are not only dissipated to heat via 
bulk viscosity, but also that energy is radiated away via 
neutrinos. This energy dissipation will be associated with 
a viscosity coefficient, the radiative viscosity (RV). We 
will also demonstrate that this RV coefficient is larger 
than the one for bulk viscosity. 

We start by studying the effect of density oscillations 
on the composition of the system, mainly defined by the 
electron fraction. From that we will find that since the 
electron fraction oscillates out-of-phase with the baryon 
density, so will the pressure, and so the integrated PdV 
work will result in a finite quantity which is the amount 
of energy of the oscillation converted to heat. This energy 
dissipation is associated with the bulk viscosity. Using 
the same formalism as for the bulk viscosity the increase 
of the energy radiated by neutrinos will be related to a 



viscosity coefficient which we call the radiative viscosity 
which wc find to be closely related to the bulk viscosity 
for the type of interactions considered here for the urea 
processes. 

Description of Oscillations — Let us begin by intro- 
ducing a periodic perturbation to the baryon density 
n{t) ~ na + Re (fee'"*), which will lead to a deviation 
from /^-equilibrium characterized by = /id — Mu ~ Me 

The density oscillations drive quark matter slightly out 
of P equilibrium, regulated by the urea processes; 



u + e 



(1) 



but not out of thermal equilibrium which is restored al- 
most immediately by strong processes. The correspond- 
ing instantaneous quasi-equilibrium state can be unam- 
biguously characterized by the total baryon number den- 
sity n and the lepton fraction Xg, 



(n„ + nrf)/3, Xe = ne/n, 



(2) 



where n„ and rid are the number densities of up and down 
quarks, while rig is the number density of electrons. (In 
the case of strange quark matter, one should also add 
the strangeness fraction Xs = riJn where ris is the num- 
ber density of strange quarks |13|.) Charge neutrality 
requires 



3^. 
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0. 



(3) 



Using this constraint together with the definitions in 
Eq. one can express the number densities and, in 
fact, all thermodynamic quantities of quark matter in 
terms of n and X,,. For the number densities, for exam- 
ple, one finds 



Tip = XpU, 



^{l + Xe)n, 7id = (2-Xe)n. (4) 



These number densities can also be expressed in terms of 
the corresponding chemical potentials, — ni{^i). In /3 
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equilibrium, the three chemical potentials are related as 
follows: fid = fJ'u + IJ-e- In pulsating matter, on the other 
hand, the instantaneous departure from equilibrium is 
described by the small parameter 



(5) 



where Sfii denotes the deviation of chemical potential /ij 
from its value in /? equilibrium. The quantity 6fi can be 
conveniently expressed in terms of the variations of the 
two independent variables 6n and SX^, 



6fi = C— + BSXe 
n 



(6) 



where, as follows from the definition, the coefficient func- 
tions C and B are given by 



dUd dfiu dfJLe 

C = rid-r. riu-r. ^e- — , 

and oriu one 



B = 



dnd dnu dn^ 



(7a) 
(7b) 



When (5/i is non-zero the two urea processes, c.f. Eq. ([T}, 
have slightly different rates. To leading order in 6fj,, we 
could write 



(8) 



The net effect of having different rates for the two pro- 
cesses is a change of the electron fraction in the system: 



n—^ = ^'5^' 



(9) 



This has the tendency to restore the equilibrium value of 
Xe- Since the rate is finite, however, the weak processes 
always lag behind the density oscillations. In order to see 
this explicitly, we substitute 5fi from Eq. (|6]) into Eq. ([9]) 
and get the equation for SX^ in a closed form. 



d{5Xe) 
^ dt 



A ( C— + B5X, 



(10) 



The periodic solution to this equation can be found most 
easily by making use of complex variables. Denoting 
SXe = Re {SXefl e'"*) , we derive the following result: 



5Xefi = 



5no 



C 



(11) 



In the last equation, the lag of the weak processes is 
indicated by a non- vanishing imaginary part of 5Xf.fi. 
Such an imaginary part controls the phase shift of the 
5Xe oscillations with respect to the oscillations of density. 

Bulk viscosity — Let us now review the derivation of 
the bulk viscosity caused by weak interactions. This has 
been done for several possible phases of nuclear matter, 
see, for example, Refs % i, i, [iMl]. 



For a periodic process, the bulk viscosity C is defined 
as the coefficient in the expression for the energy-density 
dissipation averaged over one period, r = 27r/a;, 



{£diss) 



(12) 



where v is the hydrodynamic velocity associated with 
the density oscillations. By making use of the continuity 
equation, h + nSJ ■ v — 0, we derive 



,6 \ C^"^ ( Sno 
2 \ n 



(13) 



In order to solve for C, the dissipated energy on the left- 
hand side has to be calculated explicitly. For the bulk 
viscosity, the energy dissipation is taken as the PdV work 
done on the system. 



PVdt 



(14) 



where = 1/n is the specific volume. 

The pressure oscillations around the equilibrium value 
are driven by the oscillations of its two independent vari- 
ables, i.e., the quark number density and the lepton frac- 
tion. 



dP 

5P = —5n~nC5Xe 
on 



(15) 



where C is the same as in Eq. (j7al) . In the derivation we 
took into account that ni — dP/dfii and that the total 
pressure is given by the sum of the partial contributions 
of the quarks and electrons, P = J^i Piil^i)- 

Taking into account the relation (IT5|) together with 
the solution for SX^.^ in Eq. (fTTj) . the expression (fT4)) 
becomes 



, A , I [dm 
2 V n 



(\B/n) 



(16) 



By comparing this with the definition in Eq. (jl3p . we 
finally derive an explicit expression for the bulk viscosity. 



c 



AC2 



cj2 + {XB/nY 



(17) 



For more details on the formalism and calculation of the 
bulk viscosity, the reader is advised to consult Refs. [3, 
IS, EMI. One point worth mentioning here, is that 
this bulk viscosity coefficient controls the conversion of 
the mechanical energy of the oscillation into heat (recall 
that energy dissipation is given in terms of a PdV work) . 

Radiative Viscosity — Another source of dissipation, 
which has never been investigated before, is the change 
of neutrino emissivity due to the deviation from equi- 
librium. This energy dissipation will not increase the 
temperature, instead, the energy is radiated away in the 
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form of neutrinos. We will associate a viscosity coeffi- 
cient with this energy dissipation and refer to it as the 
radiative viscosity (RV) . 

Since the change of interaction rates always leads to 
an increase of the neutrino emissivity, we can expand 
the off-equilibrium emissivity in terms of 5^^ . To lowest 
order, 



E-Eo = SE = S6fi^ 



(18) 



where Eq is the neutrino emissivity at equilibrium as cal- 
culated in Refs. [il,[i3l- Using Eqs. ^ and © we find, 
using Sfi — Re (i5/ioe*"*), that the volume oscillations are 
accompanied by an out-of phase oscillation of (5/i. 

The time-averaged change of neutrino emissivity is, 
then, 



{SE) 



(19) 



Tp{S^^) = 6 



{2TifEdEuE^Ee 



This change of emissivity takes energy away from the 
oscillation and increases the rate of energy loss via neu- 
trinos i.e. (Sdiss) — —{SE). Since this energy loss is very 
similar to the energy loss due to the bulk viscosity, it 
is justifiable to define the RV coefficient (7^) associated 
with this dissipation by using the same equation for the 
bulk viscosity c.f. Eq. (|T3l) : 



n = 2{5E)[^-^] ^-^ 



(20) 



The dimensionless ratio S jX can be found by examining 
the calculations for both the emissivity and the interac- 
tion rate. Now let us turn to the calculation of A defined 
by Eq. (|8]). Following the original approach of Iwamoto 
I16ll, we get the rate for fi decay in the following form: 



MYS" {Pd -Pu-Pe- Pp) / {Ed - ^id) [1 - / [En ~ flu)] [l-f (Ee " fie)] ■ (21) 



r 



Here, Pi and pi are the 4- and 3-momenta of the zth 
particle, respectively, and f{E) = l/(e^/-^ + 1) is the 
Fermi distribution function. The scattering amplitude 
squared is given by 16 1 



\M\ 



64G% CDs'" OciPd ■ P,){Pu ■ Pe) 



2V 



Gl cos^ 0cEuEdEpEe{l - cos ^^p) (22) 



After substituting this approximate form of |Afp, all an- 
gular integrals in Eq. ([?T|) can be done exactly. Then, 
using the dimensionless variables Xi — (Ei — fJ.i)/T (note 
that /ip = 0), we obtain 

^ —^G^P COS^ 9c fidfiufJ-eT^ / dXpSp J(xp - f), 

(23) 

where ^ = 5fi/T, and 



J{x) 



3 

n 



dxjfixj) 



S{xi+X2+X3—x) = 



By noting that Tjp{£_) — ry(— ^), we find 
17 



A = 



IStt- 



rG^cos 9casfidfJ.ufieT 



2(1 + 6=^)' 

(24) 
(25) 



In order to find the value of S we need to evaluate the 
off-equilibrium neutrino emissivity E = E[, + E^, where 
Ep and E,^ are the energy loss rates via anti-neutrinos 



and neutrinos, respectively. The integral for the energy 
loss rate due to anti-neutrinos Ep is the same as the one 
in Eq. (PT|) except that the integrand is multiplied by 
the energy of the anti- neutrinos pp, and so, the angular 
integral and the integral over the momenta of quarks and 
electrons are the same; 



MO 



^GjpCos^OcfJ.dfJ-ufieT'' I dxi^x$J{xp - ^) 



Noting that i?p(^) — Ei,{~^), we derive 



(26) 



E ^ ^Glcos''{e,)fidtiutieT'^ 



45 77r6 ITtt* fSfi 



2520 40 V T 



(27) 

The first term in this equation is the equilibrium emis- 
sivity. We finally find that 



5 3 
A ~ 2' 



(28) 



We stress that this ratio is generic to all urea pro- 
cesses, both in nuclear matter and quark matter, see, for 
instance Refs. H, 19| . 

Discussion — In this letter we discussed the effect of 
urea processes on the dissipation of density perturbations 
in the star. It was shown that they do not only contribute 
by converting energy into heat "via bulk viscosity" , but 
also that they dissipate the energy by converting it into 
an increase of the neutrino emissivity. 
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We have associated this energy dissipation with a vis- 
cosity coefficient, the radiative viscosity (RV), which, ex- 
cept for not converting the energy of the perturbation 
into heat, behaves exactly hke a bulk viscosity coeffi- 
cient: the mechanical energy of the perturbation is ra- 
diated away from the star in the form of neutrinos. We 
have shown that in the case of non-strange quark mat- 
ter, this effect is 1.5 times larger than the one of the 
bulk viscosity. This result also holds for the case of a 
proton-neutron-electron system. 

In the case of strange quark matter, the radiative dis- 
sipation is only governed by the urea processes of both 
the strange quark and the down quark. The non-leptonic 
weak interactions will clearly not contribute to the radia- 
tive viscosity (see Ref. for a detailed analysis of the 
contribution of both the non-leptonic and urea processes 
to the bulk viscosity) which may complicate the situa- 
tion since the non-leptonic weak interactions dominate 
the bulk viscosity energy dissipation for a wide range of 
temperatures and oscillation frequencies 

The radiative viscosity will contribute to the dissipa- 
tion of density oscillations in the star, these oscillations 
may arise during phase transitions as in Ref. [20|] , or due 
to giant flares in magnetars (highly magnetized neutron 
stars), see Ref. [2l|, or instabilities that result from the 
emission of gravitational waves [j^ - Eij l. The so-called r- 
mode (or rotation-dominated) instabilities might be the 
most important ones. They can develop at a relatively 
low angular velocity |27h3(J |. and therefore may be rele- 
vant for a large number of compact stars. The presence of 
viscosities relieves this situation by damping these insta- 
bilities, prohibiting the star from emitting gravitational 
waves and loosing angular momentum. 

We note that the dissipation of radial perturbations via 
an increase of the energy loss to radiation is a general phe- 
nomenon applicable to all physical processes that involve 
transfer of energy to radiation. Such processes are rele- 
vant, for instance, for variable stars, which change their 
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luminosity due to density oscillations, see Refs. |3ll. |32|. 
Examples for such stars are pulsating white dwarfes or 
asymptotic giant branch (AGB) stars. This change of lu- 
minosity may be associated with a radiative viscosity co- 
efficient responsible for dissipating these oscillations. In 
this case the energy is radiated away by photons from the 
surface rather than neutrinos. One has to be cautious, 
though, that the formalism used in this letter assumes 
averaging over one period of oscillation, which may give 
zero in the case of variable stars, so that the formalism 
might have to be changed to account for this fact. 
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